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Abstract 

In this paper we study the time evolution of the decay process for a particle 
confined initially in a finite region of space, extending our analysis given 
recently [Phys. Rev. Lett. 74 (1995), 337]. For this purpose, we solve 
exactly the time-dependent Schrodinger equation for a finite-range potential. 
We calculate and compare two quantities: i) the survival probability S(t), i.e., 
the probability that the particle is in the initial state after a time t; and ii) 
the nonescape probability P(t), i.e., the probability that the particle remains 
confined inside the potential region after a time t. We analyze in detail the 
resonant and nonresonant decay. In the former case, after a very short time, 
S(t) and P(t) decay exponentially, but for very long times they decay as a 
power law albeit with different exponents. For the nonresonant case we obtain 
that both quantities differ initially. However, independently of the resonant 
and nonresonant character of the initial state we always find a transition to 
the ground state of the system which indicates a process of 'loss of memory' 
in the decay. 
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I. INTRODUCTION 



In a recent publication JI] , we study the time evolution of the decay process in Quantum 
Mechanics for a particle initially confined in a finite region of space. We were able to solve 
exactly the time-dependent Schrodinger equation for a finite range central potential V(r) 
by considering the full resonant spectra of the system using a novel representation of the 
time-dependent Green function in terms of resonant states. 

Once we obtain the wave function -0(r, t) for our problem we calculate and compare two 
quantities: i) the survival probability S(t), that is, the probability that the particle is in the 
initial state after a time t; and ii) the nonescape probability P{t), defined as the probability 
that the particle remains confined inside the potential region after a time t. Of the two 
notions the survival probability has been more widely used, presumably because it can be 
manipulated in a more abstract way, without referring to an initially confining region of 
space for the decaying particle. One finds that the notion of survival probability has been 
commonly used in studies of quantum decay [0 and also in dynamical studies of quantum 
chaos ||. Actually the survival and nonescape probabilities correspond to two different 
notions [|]||. For instance, a particle can leave the initial state ip(r, 0) while remaining 
inside the potential. In such a case S(t) varies with time while P(t) is always unity. 

In this paper we will analyze in detail resonant and nonresonant decay. The former refers 
to the situation where the energy of the initial state is close to one of the resonances of the 
system. This case was considered in Ref. [p]] and here we extend the discussion also to the 
nonresonant decay, that is the case where the energy of the initial state is arbitrary. 

A comparison between the survival and nonescape probabilities might be of particular 
interest in connection with the problem of the nonexponential contributions to the time 
evolution of quantum decay @J^J31; in threshold effects in photoionzation in atomic physics 
H and in radioactive decay of nuclei PJTOl ; in one dimensional quantum tunneling fTT| 



m 



tunneling dynamics of squeezed states |12|], and in studies of the time evolution of chaotic 
quantum systems PJlI 
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A novel feature of our approach is that the exponential and nonexponential contributions 
to the time evolution of decay are presented in terms of the full set of complex poles and 
resonant states associated with the outgoing Green function of the system. The description 
covers the full time interval from < t < oo and in particular it allows to study the effect 
of the far away resonances on the time evolution properties of the system. 

We shall discuss the above notions for a symmetrical potential V(r) as shown in Fig. [l], 
and for simplicity, we shall restrict ourselves to s-waves. A first question is whether S(t) 
and P(t) can be obtained in an analytic and explicit fashion, and the answer is yes, if ip(r, 0) 
is available in an explicit analytic form. This discussion will be carried out for arbitrary 
V(r) and ip(r, 0) in section II using an approach by Garcia-Calderon ||. First we provide 
the stationary Green function of the problem as an expansion involving the resonant states 
and complex poles of the problem and then, with the help of the inverse Laplace transform, 
get the time-dependent Green function as an infinite sum of terms. Each term consists of 
coefficients involving the overlap of the initial state and the corresponding resonant state, 
multiplied by a function (introduced long ago by Moshinsky |TJ|]) that depends on the 
complex energy and the time. 

Once the probabilities of survival and nonescape are explicitly determined comes the 
interesting question of discussing their behavior and then comparing S(t) and P(t) for 
different times. For the resonant case we obtain that, after a very short time, both quantities 
decay exponentially and coincide with each other; however for the nonresonant case in 
general differ considerably. Independently of the resonant or nonresonant character of the 
initial state we always find a crossover to the ground state of the system which indicates a 
process of 'loss of memory'. We will also derive with more detail the main result of Ref. [|TJ, 
namely, that for long times, S(t) and P(t) behave in a different way: P(t) ~ t^ 1 while S(t) 
obeys the well known power law decay ~ t~ 3 . 

To illustrate the approach it is convenient to consider a specific example and thus in 
section III we take the case of the delta- function potential, V(r) = b5(r — R), and a simple 
analytical expression for the initial state ip(r, 0) = (3{k) sin nr. We shall in particular analyze 
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the case where the strength of the 5-function potential is very large, i.e. in our units b >> 1, 
because then one can derive an approximate analytical expression for the poles of the Green 
function. This will allow us to discuss in section IV what happens to S(t) and P(t) when 
the initial state is near resonance and compare it with the situation when the initial state is 
off-resonance. Another aspect of interest is the behavior of S(t) and P(t) at different times 
and in particular the change from exponential to a power law behaviour which we shall also 
illustrate numerically. 



II. DETERMINATION OF THE SURVIVAL AND NONESCAPE 

PROBABILITIES 

The situation that we may confront then is that of a potential V(r), as drawn in Fig. [I], 
that terminates at r = R. Our initial normalized wave function is restricted to the interval 
< r < R. We can then solve the time-dependent radial Schrodinger equation with the 
potential V(r) and the initial condition mentioned above at t — 0, and denote the solution 
by ip(r,t). With the help of the latter we can define the survival amplitude A{t) P]J4]J5|] as 



rR 

A(t) = / i/;*(r,0U(r,t)dr (2.1) 
Jo 

so that S(t) = \A(t)\ 2 is the probability of finding at time t the state ip(r,t) at its initial 
value ip(r, 0), that is, 
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S(t) 







R 

ip*(r, 0)^(r, t)dr 



(2.2) 



On the other hand the probability that the particle does not escape from the potential is 

rR 

p(t) = / ^*(r,tU(r,t)dr, 2.3 
Jo 

as the wave function ip(r,t) remains normalized if that was the case for its original value 
■0(r, 0). Thus the probability that the particle escapes is 1 — P(t). We shall denote P(t) as 
the nonescape probability. 
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To derive a procedure for determining S(t) and P(t), we require first to find a solution 
of the time-dependent Schrodinger equation for a potential that vanishes beyond a distance, 
i.e., V(r) = 0; r > R, (Fig. |1]), and the initial condition tp(r, 0). It is convenient first to 
determine the time dependent Green function, which we denote by g(r,r',t) and satisfies 
the equation 



Or 



—2 + V(r) 



d 

g(r,r,t) = i—g(r,r',t), 



with the initial condition 



g(r, r', 0) = 5{r — r'), 



(2.4) 



(2.5) 



where we take h = 2M = 1 and hence E = k 2 . 

We proceed then to find the Laplace transform of this Green function which we denote 
by g(r, r', s) given by 

("OO 

g(r,r';s) = / g(r,r' ,t)e~ st dt. (2.6) 
Instead of the complex variable s it is convenient to introduce the variable k obeying, 

s = (-^ 2 ), k = (ts) 1/2 , (2.7) 
and to define the outgoing Green function G + (r, r'; k), which is used normally in the litera- 



ture 15 through, 



G + (r, r'\ k) = ig(r, r'\ s). 



(2.8) 



^2 + k2 - V(r) 



Applying then the transform to both sides of Eq. ( p.4| ) we obtain, using Eq. ( |2.8| ), 

G + (r,r';k) = 8{r-r'). (2.9) 

Once the Laplace transform (|2.6|) is determined, we can get g(r,r',t) by the inverse Laplace 
transform over the Bromwich contour but, in terms of the variable k, the integral is over the 
hyperbolic contour Cq of Fig. |2] and we obtain 

(2.10) 

Thus our first objective will be to obtain G + (r,r'; k). 



g ( rj r ', t) = — G + (r, r'; k)e~ lk 1 2kdk. 
2tt Jc 



A. Determination of G + (r, r';k) 



Turning our attention back to Eq. fl2.4j ) we note that for r 7^ r', G + (r, r';k) satisfies 
an homogeneous ordinary differential equation, namely the Schrodinger equation of the 
problem. It is well known that G + (r, r'; k) may be written as fI5(| , 

^(r/i^- ^yH (2.1D 

where r< indicates the smaller of r and r' and r>, the larger. The function <j)(k, r) stands for 
the regular solution of the Schrodinger equation of the problem with boundary conditions 
at the origin, 



4>(k, o) = o, 



d<f)(k, r)/dr 



r=0 



(2.12) 



and the function f(k,r) is irregular at the origin and is defined by the condition that for 
r > R (since the potential vanishes beyond the distance R), it behaves as exp(ikr), namely 
an outgoing wave, which implies that for at r = R it satisfies, 



f(k,R) = e 



ikR 



df(k, r)/dr 



= ike ikR , (2.13) 

r=R 



which also defines it completely. The function W(k) represents the Wronskian, 

wik) = mr fjm. mr fjm. (2 . 14) 



It is well known [[L5[] that for a finite range interaction the outgoing Green function 
G + (r,r'; k), as a function of k, can be extended analytically to the whole complex fc-plane 
where it has an infinite number of poles, that arise from the zeros of the Wronskian Q2.14|) , 
and are distributed in a well known manner [p]5|,|17| . Purely imaginary poles on the upper 
half fc-plane correspond to the bound states of the problem, whereas purely imaginary poles 
seated on the lower half /c-plane correspond to antibound states. On the other hand complex 
poles are only found on the lower half fc-plane. To each complex pole at k n = a n — ib n , with 
a n , b n > 0, there corresponds, due to time reversal invariance, a complex pole fc_ n situated 
symmetrically with respect to the imaginary axis, i.e. = — k* n . In the examples that we 
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shall discuss later, we shall assume potentials with no attractive part so they have no bound 
states, and thus all the poles of G + (r, r'; k) will be in the lower half of the complex fc-plane. 
Hence our analysis is restricted to resonant states. 

As shown below the complex poles of the outgoing Green function are related to the 
resonant states of the system. Resonant states may be defined as solutions of the Schrodinger 



equation of the problem [|18| , |19| , |16 



d 2 u n {r) 
dr 2 



+ 



k 2 n - V{r) u n {r) = 0, 



(2.15) 



that obey the usual boundary condition at the origin, 



«n(0) = 0, 



(2.16) 



and at distances, r > R, describe a situation where there are no incident particles, namely, 
u n {r) oc exp(ik n r), and hence satisfy the outgoing boundary condition, 

du n (r) 



dr 



ik n u n (R). 



(2.17) 



r=R- 



Gamow |TB| showed that the outgoing boundary condition Q2.17p implies that the energy 



eigenvalues of the problem are necesarily complex, namely, k\ = E n 



iT n /2, where 



e n stands for the position of the resonance and T n refers to the corresponding width. The 
above boundary conditions also provide the bound and antibound solutions corresponding, 
respectively, to imaginary positive, k n = i^ m) and imaginary negative, k n = —i5 n , values of 
k n , where both 7 n and 8 n are real. The above complex eigenvalues correspond precisely to 
the complex poles of the outgoing Green function of the problem G + (r,r'; k). Moreover, it 
turns out that resonant states may also be defined as the residues at the complex poles k n 
of G + (r, r'\ k) ||2T)fl . Here we follow Garcia- Calderon to write, 

u n {r)u n {r') 



R{ r y,k n )= lim \(k-k n )G + (r,r',k)\ 



(2.18) 



provided the resonant states normalized according to 

rR 
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[ R u 2 n (r)dr + ^-K(R) = 1. 
Jo 2k„ 



(2.19) 



Let us now consider the integral || 

J_ I G + (r, r'\ z) ^ _ 

Z7TI JO Z — K 

which is taken over the contour C of Fig. ^ consisting of one large cycle, whose radius will 
eventually go to oo, and small circles surrounding all the poles of G + (r, r'; z) in the z plane 
as well as the point where z = k. As the integrand is analytic inside this contour the integral 
vanishes, as indicated on the right hand side of fl2.20|) . 

For the large circle, when \z\ — > oo, we can, in the upper half J + of the z plane, essentially 
disregard the potential V(r) in Eq. ( |2.9|) , and thus G + (r,r'; z), obtained from fl2.11|) , when 
we replace k by z, becomes 

^ je i2(r+r,) - e iz{r '- r) J, if < r < r' < R, (2.21a) 

^- 1 |e i2(r+r ' ) - e iz{r - r ' ] |, if < r' < r < R, (2.21b) 

which clearly vanishes, at least as when z = x + iy, with y > 0. For the lower half 

J_ of the z plane, and the real axis, the analysis is more complicated. Garcia-Calderon and 
Berrondo |2l| have shown though, using appropriate forms of the Born approximation, that 
G + (r,r';z) vanishes there exponentially provided (r, r') < R. Thus there is no contribution 
from the large circle in Fig. |3|, and considering the contributions from all the small circles, 
including the one around z = k, one obtains Pj22|-p5|| , 

G + (r,r'-k)= £ If^y (r,r')<R. (2.22) 

n =-oo ^n\fy r^n) 



Substitution of Eq. ( [2.22|) into Eq. (|2]9|) yields after straightforward algebra the relations 

IB 

£ Un{r ] Un{r,) =0; (r,r')<R, (2.23) 



and 



oo 

- J2 u n {r)u n {r')=5{r-r')- (r, r') < R. (2.24) 



Eqs. Q2.24Q have been also derived following an expansion of G + (r,r f ; k) in terms of inverse 
powers of k p2p3| . Notice that 

1 r 1 1i 

(2.25) 



k k n k n 



^kfi(^k kyij 2 A? 
so then using Eq. ( |2.23| ) one may write Eq. ( |2.22| ) as 



G + {r,r,k) = — ^ ( k -k n ) ' ( r ' r ) <i? - 



(2.26) 



B. Determination of g(r, r',t) 

One way to determine g(r,r',t) is deforming the contour Co of Fig. ^| to exploit the 
analytical properties of G + (r,r f ; k) using the theorem of residues. In this form one may 
obtain expansions involving a sum of exponentially decaying terms plus an integral contri- 
bution along certain path of the complex fc-plane []16|,p6| . If the potential vanishes beyond a 
distance and the notions of interest are defined within the potential region, as happens for 
the survival and nonescape probabilities defined by Eqs. ( |2.2| ) and ( |2.3| ), we can proceed as 
follows. First we deform the contour Cq to the contour C as shown in Fig. |2|. This can be 
easily done because exp(— ik 2 t) converges as k increases in the upper left quadrant of the 
A;-plane. Since we have assumed absence of bound states we then may write Eq. ( p . 1 0|) as, 

j roo „ 

g (r, r', t) = — G + (r, r'; k)e~ lk * 2kdk. (2.27) 

Z7l J-oo 

We can then substitute Eq. ( |2.26|) into Eq. (|2.27|) to obtain an expansion of g(r,r';t) in 
terms of the resonant states of the problem. The only integral we need to determine is 

2 roo q ik^t 

7T- / T—rdk } (2.28) 

Moshinsky [^] has discussed the above type of integral. It appears in the description of 
transient effects in time dependent problems in Quantum Mechanics an< ^ ^ ma y ^ e 

denoted as, 



9 



i r°° e * fc2 * 1 > 



M(k n , t) = — - -dk = -e u erfc(u), (2.29) 

where 

u = -exp(-ivr/4)A; n t 1/2 . (2.30) 

Thus the time dependent Green function g(r, r', t) takes the form first derived by Garcfa- 
Calderon ||, 

oo 

g(r,r',t)= £ Mn (r) Mn (r')M(fc n , t) , (r,r') < R (2.31) 

n=— oo 

Having obtained our basic result, we proceed to the determination of ip{r,t) for an 
arbitrary initial condition ip(r, 0). 

C. Expansion of ^>(r, t) 

If our initial state is ip(r, 0) instead of 5(r — r'), it is clear that at time t ip(r, t) is given 
by the integral 

r R 

if)(r,t)= il>(r',0)g(r,r',t)dr'. (2.32) 



o 

We define now C n and C n by the expressions 

C n = j R 'iP(r,0)u n (r)dr; C n = [* f(r,0K(r) dr, (2.33) 
Jo Jo 

where, respectively, u n {r) is normalized according to Eq. ( 2.19 ) and ip(r,0) as, 



(r,0)\ 2 dr=l. (2.34) 

o 

Finally we can write ip(r,t) as 

oo 

*l>{r,t)= £ C n u n {r)M{k n ,t); (r < R). (2.35) 

n=— oo 

The coefficients ( |2.33| ) obey some useful relations. Multiply both Eqs. ( |2.23| ) and fl2.24| ), 
respectively by if)(r, 0) and ip*(r', 0) and integrate the result from the origin up to the radius 
R, to obtain ||, 
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oo /~i /~i 



(2.36) 



and 



1 oo 

2 E = 1- 



(2.37) 



D. Expansion of the survival and nonescape probabilities 

From flOD , O and flOSD we get that 



A(*)= 2 C n C n M(k n ,t), (2.38) 

n=— oo 

It is worth mentioning that as long ago as 1951, Moshinsky had given an exact expression of 



A{t) for the one-level case, i.e., when only k±i ^ 1, see Eq. (27b) of ref. |29[. Substitution 



of Eq. ( |2.38|) into Eq. (|2.2|) leads to the following expression for the survival probability, 



s ( t ) = E E C n C n C* i C* t M{k n ,t)M*{k e ,t). (2.39) 



n=— oo £=— oo 



Using ^(r, t), as given by ( |2.35| ), one may obtain an expression for the nonescape probability 
P(t) with the help of ( p.3|) . We only need to add the definition 

I nl = [ R u}(r)u n (r)dr (2.40) 



o 



to obtain 

p ( t )= E E C n C* e I ne M(k n ,t)M*(k £ ,t). (2.41) 



oo oo 



71= — OO l= — O0 



Thus, we have the general expressions for the probabilities of survival and nonescape 
in the form of a double sum involving the products M(k n ,t)M*(ke,t) of functions defined 
by ( p.29| ), with coefficients of the type C n ,C n ,I n £ given respectively by fl2.33j ) and ( [2.4Q|) . 



Notice that the only difference between both expansions is that in P(t), I n e replaces the 
product 0*0% that appears in S(t). 
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E. Exponential and long time behaviour 



Eqs, ( |2.38| ), ( |2.39| ) and ( |2.4H ) are given in terms of M functions and therefore their ex- 
ponential behaviour is not exhibited explicitly. This can be achieved by using the symmetry 
relations between the poles on the third and fourth quadrants, k_ p = —k* to write the sums 
only over the poles k p located on the fourth quadrant and make use of the relation H6|,|28||, 



M(k p , t) = e- ik i l - M{-k p , t), (2.42) 



where M(—k p ,t) is defined as in Eq. (|2.29|) but with the argument u = exp(— in / A)k n t l l 2 



instead of fl2.30|) . Using the above relation the survival amplitude ( 2.38 ) may be written as 



A(t) = Y,C P C p e~^ t -I(t), (2.43) 
P =i 

that displays explicitly the exponentially decaying behaviour recalling that k^ = e p — iT p /2. 
The term I(t) on the right-hand side of the above equation stands for the nonexponential 
contribution and it reads, 

oo 

i(t) = E [c P c p M(-k p ,t)\ - c;c;M(-k;,t)} , (2.44) 
P =i 

where we have used that C- P C- P =C*C* and k^ p = —k*. The argument of the function 
M(-k*,t) is u = exp(-«7r/4)A;*t 1 /2_ Notice that the sums in Eqs. (g]43|) and (^441) run 
over the same set of poles. The nonexponential contribution given by Eq. ( [2.44|) appears 
as an infinite sum of terms on the same footing that the description of the exponentially 
decaying terms. Is is worth mentioning that Khalfin, who was the first author to discuss 
nonexponential contributions to decay, referred to models involving a single resonance term 
f32|| . In other treatments, the nonexponential contribution appears as a complicated integral 
term, and in general has been treated in an approximate way P J26|J2^J33[1 



If one ignores the nonexponential contribution I(t) in Eq. ( |2.43| ) and furthermore assumes 



that the initial state is very close to the sharpest resonant state, say the sth, then from ( |2.36| ) 



and ( |2.37| ) it follows that C S C S ~ 1, all other coefficients being very small, then ( |2.43| ) may 
be written as A(t) = e~ lk ^ and the survival probability becomes the well known expression, 
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S(t) 



-T s t 



(2.45) 



A similar analysis can be done for the nonescape probability. Noticing that I s 
follows from inspection of fl2.19|) and (|2.40|) , one obtains 



Pit) 



-T.t 



1, as 



(2.46) 



The equivalence of S(t) and P(t) as shown above, along the exponentially decaying region, 
has probably led to confusion regarding the notions of survival and nonescape probabilities. 

Let us now turn to the analysis of the long time behavior, as the discussion for very 
short times has been given elsewhere |28] . From Eq. ( |2.42| ) one sees that asymptotically the 
relevant terms are of the type M(—k n ,t) where k n stands for either k p or k*. Indeed from 
the definition of the function M, given by (|2.29|) , one may get the asymptotic expansion of 



exp(-u 2 )erfc(w) to obtain 
M(-k n: t) 



+ 



(2.47) 



2(7Ti) 1 / 2 VM 1/2 / 4(tu) 1 /2 \klt 3 /\ 

Substitution of Eq. ( |2.47| ) into (|2.44| ) allows to write I(t) as a sum over terms that go like 
inverse powers of time. One may see, however, that the coefficient proportional to t^ 1 ^ 2 , 
is identical to Eq. ( [2.36|) and therefore it cancels out exactly. Consequently the survival 
amplitude ( |2.43| ) may be written as, 



A(t) « E 



CpCpe-^ 



4(7Ti)V2 



CpCp 



CpCp 

U3 



£ 3/2 



(2.48) 



The above equation exhibits the crossover from exponential to a power law behaviour. Using 
( p.48| ) one sees that the survival probability (|2.39| ) behaves at long times as S(t) ~ t~ 3 . 

Let us now consider the long time behaviour of the nonescape probability. One sees from 
Eq. (|2.31|) into Eq. ( [2. 3D for P(t), that the resonant expansion of g(r, r', t) is coupled through 
the integration over r, with that of g*(r,r',t). This originated the integrals I n £ defined by 
( p.40| ). Hence, when the asymptotic expansion of the M functions ( |2.47|) is introduced into 
( p.41| ), the leading contribution is proportional to t" 1 and includes terms of the type 
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OO OO / (~t*(-1 T \ -1 
r S \ rh r fh S J i> 

that are different from ( |2.36| ) and hence do not cancel. In other words, asymptotically, 
Pit) ~ £ -1 . One concludes therefore, that the survival and the nonescape probabilities 
behave with a different inverse power of time for long times Q. 

The previous discussion refers only to resonant states, however it can be easily extended 
to include bound states 0. In such a case the long time behaviour of the survival and 
nonescape probabilities has a nonvanishing asymptotic value due to the oscillating character 
of the time dependence of the bound states. For instance, if in addition to resonant states 
there is a bound state, as time goes to infinity the survival probability reads, 

S(t) = \C b \ 2 , (2.50) 



where C& and C& = C£ stand for the overlap intregrals defined in Eq. (|2.33|) with the u n (r) 



replaced by the bound state function tt&(r). The above nonvanishing asymptotic value for 
the survival probability is sometimes referred to in the literature as population trapping || . 
Proceeding in a similar fashion as above for the nonescape probability yields, in the limit 
that the time goes to infinity, that 

Pit) = \C b \ 2 . (2.51) 

One sees that both S(t) and P(t) provide at very long times a similar result when bound 
states enter into the problem. 

III. EXAMPLE 

We shall now consider a problem that illustrates numerically the behaviour of S(t) and 
P(t). We take V(r) delta-function potential 

V(r) = b5(r- 1), (3.1) 

with strength b and where we choose R = 1. For the initial condition we take 
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</>(r,0) 



(3{k) sin/tr, if < r < 1, 
if r > 1, 

where k is an arbitrary real parameter and, for normalization, we require that 



(3.2) 



(3{k) = V2[l - (2k)- 1 sin 2k]- 1 . 



(3.3) 



V(r,0) 



The form (|3.2|) was suggested by the fact that if b — > oo, the normalized stationary states 
for the problem are 

\/2 sin mnr, if < r < 1; m = 1, 2, 3, . . . 

(3.4) 

ifr>l, 

and so an interesting question could be what happens at any time t if at t = we consider 
the above initial state. This state implies that initially we are very close to a resonance, 
and since we would like to consider also the case of an initial state situated between two 
resonances, we use the more general form ( |3.2|) . 



From Q2.15|) we immediately see that 



u n (r) = d n sink n r, 



(3.5) 



where d n and the equation satisfied by k n will be given below. Using (|3.5| ) we see that Eq. 
( p.l?| ) has to be modified to take into account the fact that the delta-function potential 
is discontinuous at r = R = 1. Hence for this example instead of Eq. ( |2.17| ) we have the 
condition 

ik n u n {\) - l^j^-} = &u»(l). (3.6) 

Substitution of Eq. (|3.5| ) into the above expression yields the equation for the complex 
eigenvalues of the problem, namely, 



2ik n + b(e 2lkn - 1) = 0. 



(3.7) 



The solutions to the above trascendental equation correspond also to the poles k n , (n = 
±1, ±2, ±3, . . .) of G + (r,r f ; k) and in this case all of them are in the lower half of the k 
plane with k_ n = — fc*. 
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The coefficient d n in Eq. (|3.5| ) may be obtained from the normalization condition given 
by Eq. ( |2.19| ). Using ( |3.7| ) it may be written as, 

_ r 2(b-2 ikn ) i 1/2 

n ~ \{\ + b-2ik n )\ • (3 - 8) 
In order to calculate the survival and nonescape probabilities, given respectively by ( |2.39|) 
and ( p.41| ), we need to determine also the coefficients C n and I n g using (|3.5| ) for u n (r) and 
( p.2|) for ip(r, 0). Since the initial condition ip(r, 0) = /3(fi;)sinKr depends on the parameter 
k, we shall denote C n as C n (K). Hence from ( |2.33|) we obtain 



q ihn 

C n (K) = {3(K)d n -r-r\ r^T- [k COS K - i(k n + lb) sin k] . 

Note that as our initial state ( |3.5| ) is real C n = C n . On the other hand from (2.4U) and 



(3.9) 



d n d}k n k} j(k*-k n ) 
e , 



6 2 (/c| k n ) 

where for all of these coefficients we have made use of the fact that, from 

2ik n _ 2k n + ib 
ib ■ 



(3.10) 
, k n satisfies 
(3.11) 



Substitution of the coefficients given by (|3.8|), (|3.9|), and fl3.10| ) in ( p.39| ) and ( p.41| ), leads 
to explicit analytic expressions for the probabilities of survival S(t) and nonescape P(t) that 
we now proceed to discuss. 

An interesting case refers to sharp isolated resonances. For not very highly excited 
resonances, this case can be achieved by taking a large value for the intensity of the potential, 



namely, in our units, b ^> 1. In that shown in ref. J28|, the poles of G + (r,r'; k) are 

given approximately by 



1 



nir\- 



(3.12) 



A. Comparison of S(t) and P(t) 

We shall analyze the survival and nonescape probabilities when b 3> 1 and k = run. The 
complex energy E n associated with k n of (|3.12|) is approximately given by 
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3^3 



.n it 



E n = k 2 n ~n^-2i—, (3.13) 

so that the separation of the resonances, i.e. the real part of the difference E n+ % — E n , 
is much larger than the width of one of them, i.e. the imaginary part of E n . The choice 
k = mir implies that the initial state ( |3.2|) is associated with an energy k 2 = m 2 ir 2 and thus 
it is very close in energy to the corresponding resonance as follows from (|3.13| ) with n = m. 
Since /3(m.7r) = (2) 1//2 , we see that the formula (|3.9| ) for C n (nm) becomes 



2 7^777. 

C n {mn) » — — — ; n ^ m, (3.14) 

0(777/ — n z ) 

that has a large factor b in the denominator. On the other hand for n = m, we have from 
( gTTg ), that (m 2 7T 2 - k 2 m ) « {2m 2 -n 2 /b) and thus 

C m (m7r) « 1. (3.15) 



We can also analyze for b ^> 1, and 6 much larger than n and £, the behaviour of J n ^. 
We get from ( ^T0|) and fyj) that 

2in£ir 



b 2 (£~n) 



n ^ 



(3.16) 



and 



7 « 1 



(3.17) 



Thus, if 6 ^ 1, and also 6 much larger than m, and the initial state is very close to the 
rath resonant state, we see from fl3.14| ), ( |3.15| ) and ( |3.17 ) that the survival and nonescape 
probabilities essentially reduce to 



S(t) « P(t) « e 



-T m £ — t/r 



(3.18) 



with r defined as 



4m% 3 
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(3.19) 



Hence when we are close to a narrow resonance, the survival and nonescape probabilities 
coincide for a number of lifetimes until the contribution of the resonance with the largest 
lifetime dominates the time evolution of S(t) and P(t), as will be shown in the numerical 
examples presented in the next section. 

On the other hand when the initial state is between resonances, i.e. k = (m + l/2)n, 
inspection of Eq. (|3.9|) yields 

C n ((m + 1/2)tt) « -f— n ± m, (3.20) 

(m 2 — n 2 )7T 

and for n = m 

C n ((n + l/2)7r)«--, (3.21) 

7T 

which indicates that there is no a single coefficient C n that dominates over the expansions 
for S(t) or P(t). Hence there is no reason to expect initially an exponential decay law 
behaviour. In the next section we also exemplify this situation. 



B. Numerical results 

In order to analyze S(t) and P(t) we need to specify two parameters: the strength b 
of the delta-function potential and the wavenumber k for the initial condition. In what 
follows, we will take the fixed value b = 100. On the other hand, we have to find the poles 
k n in the complex fc-plane, that satisfy Eq. (|3.7|) . In order to do so, we use the well-known 
Newton-Raphson method to locate the poles in the fourth quadrant k n (n = 1, 2, . . . , N); the 
remaining poles lie on the third quadrant and are given by /c_ n = — k* n . We have considered 
in our calculations 1000 poles (N=500); those on the fourth quadrant are depicted in Fig. 
[|. Notice that the real part of k n is much larger than the imaginary part of k n and, as a 
consequence, we are in a regime in which the real part of the complex energy E n is much 
larger than its imaginary part, that is, we are dealing with sharp isolated resonances. The 
poles in Fig. |3] depend on the parameter b, and when b » 1 they are given approximately 
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by Eq. Q3.12|) . In fact, we use ( |3.12|) as an input for the Newton- Raphson algorithm, in order 



to obtain the exact location of the poles. When b — > oo, the poles "move" to the real axis in 
the k— plane and are given by k n = mr; this limiting case corresponds to the bound states 
of a particle in a box. 

Once we have obtained the analytical expression for the survival S(t) and nonescape 
P(t) probabilities (Eqs. (gjg) and (fTIlD), we can analyze both quantities numerically. As 



mentioned above we will take b = 100. In the remaining figures we plot the logarithm of 
S(t) as a solid line and the logarithm of P(t) as a dashed line. In Fig. 5 we show In S(t) 
and lnP(t) as a function of time, when the initial condition has a value k = 5tc. This 
value corresponds to an initial state close to the fifth resonance and hence refers to the time 
evolution of an initial state at a higher excitation energy to those considered in Ref. [1], 
where the initial states were considered, respectively, near the first and second resonances 
of the system. The fifth resonance has a lifetime t$ = I/T5, where is its width. We will 
take the time in Figs. 5 and 6 in units of r 5 . 

We see from Fig. 5a that at the beginning of the decay process, both quantities coincide 
and both decay exponentially with a lifetime r 5 , as expected from the choice of the initial 
condition. As is well known ]28| the exponential decay is not valid for very small times, 



i.e., times much smaller than a lifetime, however this deviation cannot be appreciated in 
our Figures due to the scale. In Fig. 5b we can see that after a time ~ 10r 5 , S(t) and P(t) 
separate from each other; S(t) decays exponentially with a slope T 5 for about 20 lifetimes 
(t ~ 20r 5 ). After that time, something interesting happens: S(t) starts to oscillate and then 
decay exponentially once again, but this time with a lifetime T\ which corresponds to the 
ground state of the system. The nonescape probability P(t) , after a time m IOT5, change the 
slope of the exponential decay from T 5 to r 1; but in a smooth way without any oscillations. 
Finally, in Fig. 5c, we show the behavior of S(t) and P(t) from the beginning of the decay 
process until very long times (t/rs = 6000) in order to analyze the asymptotic decay. We 
can clearly see three separate regions for both S(t) and P(t) : the exponential region with 
slope T5 at the beginning; the crossover to an exponential decay with slope Ti corresponding 
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to the ground state and; another crossover (for 77V5 « 4000) to a power law decay in which 
S(t) ~ t~ 3 and P(£) ~ t -1 . Notice that P(t) is always a smooth function and that S(t), on 
the other hand, fluctuates during the transitions from one regime to another. 

Before we analyze the nonresonant case, let us try to understand the results depicted in 
Fig. 5. First of all, when we choose an initial condition ip(r, 0) such that k = run, with m 
an integer, we are very close to a resonance given by the pole k m . In this case we have a 
leading term for S(t) and P(t) due to the resonant denominator of the form (m 2 7r 2 — fcjjj . 
The presence of this leading term guarantee an exponential decay with a lifetime r m = T" 1 , 
associated with the complex energy E m = e m — iT m /2. Furthermore, as discussed above, 
when k = mir and b » l,P(t) and S(t) coincide. This is precisely what we found in the 
example shown in Fig. 5a, when m = 5. 

Now, in order to understand the crossover to the ground state with a lifetime T\ = T^ 1 , 
we have to remember that in our example there exist an ordering t% > r 2 > . . . > tn, that 
is, the ground state is the state with the longest lifetime and therefore it dominates the 
decay process at long times. In terms of the widths, the first resonance is the thinnest. The 
fluctuations obtained for S(t) in the interval 20r 5 to 2OOT5 in Fig. 5b, correspond to the 
interplay of all the states between m = 1 and m = 5. In addition the the cases m = 1 and 
m = 2 considered in Ref. Q, we have also analyzed the case in which m — 3, (not shown 
here) where S(t) decays exponentially with a lifetime r 3 at the beginning and, during the 
crossover to the ground state, it oscillates periodically with a period given by (e 2 — ei) _1 . 

Finally, as discussed for the general case of a finite range potential (section II), the 
asymptotic behaviors of S(t) and P(t) are governed by the asymptotic form of the func- 
tion M(k n ,t) and the sum rules given by ( |2.36| ) and (|2.37Q . The result is that we obtain 
numerically for the delta-function potential that S(t) ~ t~ 3 and P(t) ~ t -1 , as expected. 

Let us now turn to the nonresonant case in which k is not a multiple of tt. We choose 
k = 5.57T, which corresponds to an initial state situated between two resonances, i.e., m = 5 
and m = 6. In Fig. 6 we show In S(t) and lnP(t) as a function of time, in units of r 5 . Once 
again, the dashed line correspond to lnP(t) and the solid line to In S(t). In Fig. 6a we can 
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see that at the beginning of the decay process neither S(t) nor P(t) decay exponentially. P(t) 
is a smooth curve but S(t) fluctuates much more strongly than the resonant case during the 
entire process. For longer times (t/r 5 100), both quantities start to decay exponentially 
with a lifetime Ti, as can be seen in Fig. 6b, even though we are in the off- resonance case 
characterized by k — 5.57T. Finally, in Fig. 6c we can see a crossover, for very long times 
(t/r 5 PS 3000), to a power law decay. Once again, in this regime S(t) ~ t~ 3 and P(t) ~ t -1 . 

Let us try to understand the differences between P(t) and S(t) from a physical stand- 
point. The nonescape probability P(t) is, by definition, the probability that a particle 
remains confined inside a certain region. Thus, the particle can change from the initial state 
to another state inside the potential and this change does not affect P(t). The only way in 
which P(t) can decrease is when the particle leaves the potential region, say by tunneling. In 
the problem we are analyzing, once outside, the particle cannot enter this region again and, 
therefore, we expect a decreasing smooth curve describing the nonescape probability. On 
the other hand, the survival probability S(t) refers to the probability that a particle remains 
after a time t in its initial state, where the latter is confined inside the potential region. 
Suppose that a particle is initially between two resonances. Before the particle leaves the 
potential region, it can make transitions from one state to another and return to the initial 
state an so on. All these transitions are manifested by the rapid oscillations of the survival 
probability. 

A result of our analysis that seems of interest to stress here is that independently on 
the conditions for the initial state eventually the resonant state with the longest lifetime 
predominates over the time evolution of decay. Apparently the system looses memory on 
whether the initial state is close to the mth or the nth resonance, or between the mth and 
the nth resonances. The last stage of the decaying process proceeds always exponentially 
before the croosover to the power law decay takes over. Indeed, at the beginning of the decay 
proccess the differences in the time evolution for different initial states may be dramatic, 
as shown for example by Figs. 5a and 6a, and also by Figs. 5b and 6b. However once the 
transition to the first resonance occurs the decay out of this state and the crossover to the 
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power law behaviour becomes similar as shown by Figs. 5c and 6c. 

The above situation may be understood in general by inspection of Eq. ( |2.37 ) in view of 



( |2.39| ) and (|2.41| ) for S(t) and P(t), respectively, and using the relation ( |2.42| ) to exhibit the 
exponential dependence of the above definitions explicitly. The initial state determines the 
values of the different coefficients C n C n and the main point is to realize that for the longest 
lived state, say the sth, even if the coefficient C S C S is small, the product C S C s exp(—ik^t) 
will eventually dominate over all other exponentially decaying terms. This follows because 
the other decaying terms, having smaller lifetimes, go to zero much faster. 

IV. CONCLUDING REMARKS 

In this work we have made use of an exact expansion for the retarded time dependent 
Green function in terms of resonant states to study the time evolution of an initially confined 
arbitrary state. Our approach allows to treat on the same footing both the exponential and 
nonexponential contributions to the time evolution of decay. We have compared the notions 
of the survival and nonescape probabilities and found that in general they exhibit a different 
behaviour with time. For the case of an initial state close to a resonance the above notions 
coincide along the exponentially decaying region. However, if the initial state is not close to 
a resonance these probabilities exhibit a nonexponential behaviour. In the above two cases, 
we have shown that eventually the resonance with the longest lifetime predominates in the 
decay process which may be seen as a 'loss of memory' of the initial condition. This occurs 
in an exponential fashion before the onset to a power law decay is established, i.e. S(t) ~ t~ 3 
and P(t) ~ t~ x . However in the case that bound states are present both quantities tend to 
a constant, i.e., S(t) ~ P(t) ~ \Cb\ 2 - 

The above formulation may be generalized to bidimensional systems, in particular to 
systems with arbitrary geometries [Q. Also, our formulation may be of interest in studies of 
the tunneling dynamics of squeezed states JT2J and in the time evolution of chaotic quantum 



systems, where nonexponential decay have also been obtained P5|-P^]. Here it could be of 
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interest, for example, to investigate the role of the power law decay for the case of quantum 
systems that behave chaotically in the classical limit. Furthermore, our results might be 
of interest in the efforts to verify experimentally the nonexponential contributions to decay 
[/J. For nuclear radioactive decay it has been reported H no deviation from exponential 
decay up to 45 half-lives for the case of Mn 56 . A more promising possibility to observe 
nonexponential contributions to decay could be in photodetachment near threshold in atomic 
physics [fjj, [38l , |3~9|l . The analysis of the present work could be of interest in that field because 
negative ions in photodetachment near threshold are usually modelled by a single electron 
in a short range potential. 

This work was partially supported by DGAPA-UNAM No. IN-101693. We thank K.W. 
McVoy for useful discussions. 
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FIGURES 

FIG. 1. Short range potential V(r) between r = and r = R. The dashed areas indicate 
schematically two resonant states, which are centered around an energy Ei and have a width Tj. 

FIG. 2. The Bromwich contour C Q in the k— plane for the integral given by eq. (2.10), which 
can be modified to the real axis C, due to the factor exp(— ik 2 t), which forces the integral to vanish 
on the portion of the circle indicated by the dashed line. The complex poles of the Green function 
(2.11) found on the lower half k— plane are also indicated. 

FIG. 3. The contour C in the z— plane for the integral given by eq. (2.17), consisting of one 
large circle and small circles surrounding all the poles of G + (r, r' , z) as well as z = k. 

FIG. 4. Location of the first 500 poles k n , given by eq. (3.7), in the complex A;— plane. Only 
the fourth quadrant is shown and notice that Re(k) » Im(fc). 

FIG. 5. Logarithm of the survival (solid line) and non-escape (dashed line) probabilities as a 
function of time, for b = 100 and the resonance case k = 5tt. The time scale is given in units of T5. 
(a) At the beginning of the decay process both quantities coincide and decay exponentially with 
a lifetime T5. The slope is given by = I/V5. (b) After a time ~ IOT5, S(t) and P(t) separate 
from each other, and both probabilities start to decay with a lifetime ti, which corresponds to the 
ground state. Notice the oscillations for S(t) during the crossover from the excited state (slope T5) 
to the ground state (slope Ti). (c) For very long times (6000 T5) we can clearly see three separate 
regions for both S(t) and P(t): the exponential region with slope Ts at the beginning; the crossover 
to the ground state with slope Ti; and finally, another crossover to a power law decay, in which 
S(t) ~ t~ 3 and P(t) ~ r 1 . 
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FIG. 6. Logarithm of the survival (solid line) and non-escape (dashed line) probabilities as a 
function of time, for b = 100 and the nonresonance case k = 5.57T. The time scale is given in units 
of T5. (a) In contrast to the resonance case (see Fig. 5a), neither S(t) nor P(t) decay exponentially, 
and S(t) displays a fluctuating character, (b) Although we are in the off-resonance case, after some 
time, S(t) and P(t) start to decay with a lifetime ti, which corresponds to the ground state, (c) 
For very long times (« 4000 T5) we can distinguish three separate regions for S(t) and P(t): a non 
exponential region at the begining; the crossover for the ground state with slope Ti; and finally, 
the transition to a power law decay, in which S(t) ~ t~ 3 and P(t) ~ t -1 . 
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